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Abstract 

> 

Gluon-induced contributions to the associated production of a Higgs and a Z boson 
are calculated with NLO accuracy in QCD. They constitute a significant contribution 
to the cross section for this process. The perturbative correction factor (if-factor) 
is calculated in the limit of infinite top-quark and vanishing bottom-quark masses. 
The qualitative similarity of the results to the well-known ones for the gluon-fusion 
process gg — >■ H allows to conclude that rescaling the LO prediction by this if -factor 
t-h leads to a reliable NLO result and realistic error estimate due to missing higher- 

order perturbative effects. We consider the total inclusive cross section as well as a 
• i-h scenario with a boosted Higgs boson, where the Higgs boson's transverse momentum 

^ is restricted to values px,H > 200 GeV. In both cases, we find large correction factors 

3 K w 2 in most of the parameter space. 
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1 Introduction 



With the recent observation of a new particle at the LHC [T|[2] and the related evidence at 
the Tevatron [3] , efforts to determine its identity are of highest priority. Among the most 
important observables are the total and differential cross sections. First measurements of 
these quantities indicate that the new particle is indeed the long-sought Higgs boson of the 
Standard Model (SM). In order to definitely confirm or exclude this hypothesis, accurate 
measurements and corresponding precision calculations of the cross section in the various 
production modes are required. 

The current theoretical knowledge of the SM cross sections is in general quite impressive 
and documented in Refs. (4|[5] . Subject of the current paper is a particular contribution 
to the so-called Higgs-strahlung process pp — > HV (V = W,Z). While it has been a 
major search mode for Higgs bosons at the Tevatron, it used to be considered of minor 
importance at the LHC due to its small cross section and large background. However, 
it belongs to the channels that were analysed by the ATLAS and the CMS experiments 
already with the first data. The signal-to-background ratio for Higgs-strahlung can be 
significantly enhanced when cutting on events where the Higgs boson is produced at large 
Pt,h [6]- 

The leading-order (LO) cross section for this process can be written as a convolution of 
the cross section for the Drell-Yan process pp — > V* with the decay rate for V* —> HV, 
where V* denotes an off-shell gauge boson of momentum k: 

a HV,DY (pp ^RV) = J dk 2 a DY (pp -» V*) dT{V * d ^ 2 RV) ■ (1) 

This relation holds exactly through next-to-leading order (NLO) QCD, i.e. 0(a s ), and 
approximately through next-to-next-to- leading order (NNLO). The QCD effects of Eq. (TTj) 
are therefore strongly dominated by the Drell-Yan corrections to the cross section a , 
they are known through NNLO QCD for the total HW/HZ cross sections [7]-[9j, and for 
HW production also differentially [To]. Typical Feynman diagrams for the Drell-Yan type 
contribution are shown in Fig.[T](a-f). They contribute to the cross section at order g 4 a™ 
(n = 0, 1, 2) and increase it by about 30% with respect to LO. Here and in what follows, 
a s = g|?/(47r), with g s the strong and g the weak coupling constant. 

Apart from the Drell-Yan-like QCD corrections at NNLO, there are top-loop-induced con- 
tributions such as the ones shown in Fig.[T](g-j). Their interference with the LO and the 
real-emission NLO amplitude is of order \ t g 3 a 2 , with A t the top Yukawa coupling, and 
their numerical impact is at the percent level 11 . 

In contrast to the NLO QCD and dominant NNLO QCD corrections, electroweak (EW) 
corrections do not respect a factorization into Drell-Yan-like production and decay, since 
irreducible (box) corrections to qq^ —> HV already contribute at NLO. The NLO EW 
corrections have been evaluated in Ref. |12| for the total HV cross sections, where they 
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amount to —(5—10)%, and in Ref. [13] for differential distributions as part of the HAWK 
Monte Carlo program, which fully includes all decays and off-shell effects of the weak boson 
V = W,Z. In distributions the EW corrections can grow to —(10—20)%. As suggested 
in Ref. [14] , NLO EW and Drell-Yan-like NNLO QCD corrections can be conveniently 
combined in factorized form, where the EW corrections modify the QCD prediction by a 
relative correction factor that is rather insensitive to the parton luminosities. 



Recently, QCD corrections to the H — > bb decay have been considered as well 15 . These 



final-state corrections should be carefully taken into account in the Higgs reconstruction. 

In this paper we focus on another type of contribution which is specific to HZ production, 
namely gluon fusion, mediated by top- and bottom-quark loops. Typical diagrams of 
this channel are shown in Fig. [2] Owing to the initial-state gluons, it cannot interfere 
with the LO amplitude and therefore contributes to the cross section at order X^g 2 a^. 
For Mu = 125 GeV, at leading, i.e., one-loop order it amounts to about 4% (6%) of the 
total Higgs-strahlung cross section at the LHC with 8TeV (14TeV) |9|. Since it has no 
lower-order correspondence, it is separately gauge invariant and IR and UV finite. The 
two initial-state gluons lead to a rather strong renormalization and factorization scale 
dependence of about 30%, thus increasing the theoretical uncertainty of the HZ relative 
to the HW process, where the gg channel does not exist (at this order). Experience from 
the gluon-fusion process gg — > H shows, however, that the LO scale uncertainty drastically 
underestimates the actual size of the higher-order corrections. Owing to the similarity of 
the gg — > H and the gg —> HZ processes in their QCD structure (same initial states and 
colour structure, both loop-induced), we expect a similar phenomenon in the latter. 

The goal of the present paper is to improve on the theory uncertainty of the gg — > HZ 
process by calculating its NLO QCD corrections. Note that they are of order and thus 
formally contribute to the N 3 LO corrections of the Higgs-strahlung process. Technically the 
described NLO calculation involves massive, multi-scale two-loop diagrams that are beyond 
present calculational techniques, so that we are forced to employ asymptotic expansions 
in the limit of a large top-quark mass. We note that the same strategy was already 
successfully applied to the calculation of NLO corrections to the related process of Higgs 
pair production via gluon fusion, gg — > HH [16] . 

The paper is organized as follows: In Section[2] we briefly outline the problem, before 
describing the details of our calculation in Section[3l Our numerical results are discussed 
in SectionHJ and our conclusions given in Section[5l 
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2 Outline of the problem 



2.1 Leading order 

At LO and in covariant gauge, the Feynman diagrams contributing to the gluon-induced 
Higgs-strahlung process can be divided into three types, shown in Fig. [2] 

(a) Box diagrams for gg — > HZ: Only massive quarks run in the loop due to the propor- 
tionality to the respective Yukawa coupling. Note, however, that these graphs tend 
to zero also in the heavy-quark limit. 

(b) Triangle diagrams for gg — > Z* — > HZ: Owing to Furry's theorem, all contributions 
from vector couplings compensate each other, so that only the axial-vector coupling 
of the Z boson needs to be taken into account. Since the axial- vector coupling is 
proportional to the third component of the weak isospin of the quark (=t^), the 
contribution of a single quark generation vanishes in the equal- mass case. Assuming 
massless quarks in the first two fermion generations, this leaves a non-vanishing 
contribution only from the third generation. The amplitudes tend to zero in the 
heavy-quark limit. 

It is interesting to note that only the longitudinal part of the Z-boson propagator 
contributes, while all contributions of the transverse part vanish. This consequence 
of the Landau- Yang theorem |l7 18 can be used at NLO to facilitate the calculation 



significantly, as will be described below. 

(c) Triangle diagrams for gg — > G° — > HZ: Only the massive-quark loops contribute 
here, where G° is the would-be Goldstone boson partner to the Z boson. The graphs 
are both proportional to the respective Yukawa coupling and to the third component 
of the weak isospin of the quark and tend to a constant in the heavy-quark limit. 

While the box diagrams (a) are gauge-parameter independent in the R^ gauge, both the 
vertices (b) and (c) depend on the gauge-parameter of the Z boson. The sum of (b) and 
(c) for each quark generation is, of course, gauge-parameter independent. 



The full result for the LO amplitudes for the process gg — > HZ can be found in Ref. [19] 
the hadronic cross section can be easily obtained using the program vhOnnlo [9}|20]. We 



have rederived the LO cross section with the full dependence on the top- and bottom-quark 
masses as a basic ingredient of our NLO calculation. 



2.2 Next-to-leading order 

The Feynman diagrams for the NLO QCD corrections to the gg — > HZ process are obtained 
from the LO gluon-fusion diagrams shown in Fig.[2](a-c) by attaching virtual and real 
gluons and quarks to internal and external quark and gluon lines in all possible ways: 
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Figure 1: Representative diagrams to hadronic HZ production of Drell-Yan type up to 
NNLO (a-f) and non-Drell-Yan-like NNLO graphs with Higgs radiation off top-quark loops; 
both types of corrections (up to NNLO) are not considered in this publication. 
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(1) H 



Figure 2: Representative diagrams to hadronic HZ production via quarkdoop-induced 
gluon fusion. It is understood that crossed diagrams have to be taken into account as 
well. 
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• For the real corrections, this results in triangle diagrams with two massive external 
momenta, Fig.[2](g,h), box diagrams with one or three massive external momenta, 
Fig.[2](i) and (d,e), and pentagon diagrams with two massive external momenta, 
Fig. [21(f) (counting off-shell gluons as massive lines). Note that Fig.|2](e) is a crossed 
version of Fig.[T](j), for example; as pointed out above, it can interfere at 0(Xtg 3 a 2 ) 
with the NLO real-emission Drell-Yan type amplitude, which has already been taken 
into account in Ref. [IT) . In the present paper, we work at 0(X^g 2 a 3 ) and need to 
evaluate the square of such terms. 

• For the virtual corrections, we encounter two- loop vertex and box diagrams, 
Fig.[2](j,l,m), as well as one-particle-reducible diagrams with two one-loop triangle 
insertions, Fig.[2](k). 



3 Details of the calculation and effective-field-theory ap- 
proach 

While the majority of the integrals could be calculated using well-known techniques, a 
general result for the massive double-box integrals shown in Fig.[2](j) is beyond current 
technology. However, motivated by observations made in Refs. 2l]-27], for example, we 



follow a strategy that has been successfully applied to higher-order corrections to Higgs 
production via gluon fusion. Instead of calculating the Feynman integrals in full generality, 
we determine the perturbative correction factor 

NLO 

in the limit of infinite top-quark and vanishing bottom-quark masses (referred to as "effec- 
tive theory" in what follows). For the gluon- fusion process, both inclusive and differential, 



it turns out that this factor is rather insensitive to the top-quark mass effects 21-27 



Using asymptotic expansion of Feynman diagrams [28, 29 , the heavy-top limit can be 



interchanged with the loop integration, which simplifies the calculation enormously. 
(i) LO amplitude 

At LO, the diagrams with top-quark loops reduce to vacuum diagrams (integrals with 
vanishing external momenta) in the effective theory. Because already at LO the loop 
integrals are UV divergent, some care is needed in the calculation of the Dirac traces that 
involve the matrix 75. Both at LO and NLO, we consistently use the 't Hooft-Veltman 
scheme 1 30 , 31 , where 75 anticommutes with the first four, but commutes with all other 
Dirac matrices. In practice, we insert 75 = — j|e AtI /p CT 7 M 7 1/ 7 p 7 cr (e 0123 = +1), keep the 
e-tensor outside of the L>-dimensional integration, and project onto four dimensions only 
after all divergent terms have cancelled among each other. 
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The result for the LO amplitude .Mo and its polarization- and colour-averaged square is 



"s« xab ^ m ^ ^p H -e z 
M z 



>Io = — 2 2u J e(ei,e 2 ,pi,P2 



a 2 a 2 



(3) 



with Mh the Higgs mass and Mz the Z mass, a and a s the electromagnetic and the 
strong coupling constants, = 1 — the sine of the weak mixing angle, and A(a, 6, c) = 
a 2 + 6 2 + c 2 — 2afr — 2ac — 26c. In Eq. ^ we make use of the momentum assignment 
g a (pi) + g b {p2) — > H(pn) + Z(pz), with a, 6 denoting colour indices and s = (pi + P2) 2 
the usual Mandelstam variable. The polarization vectors £i(pi) (i = 1,2) correspond 
to the respective incoming gluons, and £^{pz) to the outgoing Z boson. The shorthand 
e(ei, £2,^1,^2) stands for the contraction of the e-tensor with the 4- vectors in the argument. 

(ii) NLO virtual corrections 



The large- mass expansion [28, 29 , formulated in terms of the "method of regions", states 
that the loop integrand is to be Taylor-expanded in all relevant regions of loop momenta, 
and the final result of each diagram is the sum over all those expansions. If a single 
internal mass m is considered large, a loop momentum can either be large, \q^ | ~ m, 
or small, <C m, while external momenta pf are always small, \p^\ <C m. In this way, 
the virtual NLO (two-loop) diagrams with top-quark loops reduce to either massive two- 
loop vacuum integrals or products of massless one-loop triangles with massive one-loop 
vacuum integrals. As a result of the large-mass expansion, all one-particle irreducible two- 
loop graphs involving top-quark loops vanish, except for the ones with corrections to the 
ggG° vertex. Note that the latter as well as the axial vector part of the ggZ vertex each 

to restore chiral symmetry in the massless- quark 
30,31 for 75, in the following denoted as 5Z§ and 



receive an anomalous counterterm 32 
limit in the 't Hooft-Veltman scheme 
6Z§, respectively. 

Since we set m q = for q ^ t both in the propagators and in the Yukawa couplings, the set 
of diagrams to be evaluated with internal massless quarks does not contain any two-loop 
box diagrams. The only genuine two-loop integrals with non-vanishing external momen- 
tum are three-point functions as the one shown in Fig.[2](l). As their LO counterparts at 
one loop, their contribution to the cross section vanishes when working in Landau gauge, 
where the Z propagator D^ u (k) is proportional to the polarization sum f° r a physical 
vector particle Z* of mass k 2 , 

D^(k) oc -gr + ^ = E e (A)(^(A) W • (4) 

A 

The ggZ* subamplitude therefore corresponds to the decay of a massive into two massless 
vector particles which is forbidden due to the Landau- Yang theorem [17, 18 . 
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In summary, the one-particle-irreducible two-loop diagrams comprise only non-vanishing 
contributions from top-quark loops in the ggG° vertex and massless quark loops in the ggZ* 
vertex. Since the contribution of the latter vanishes in Landau gauge for the Z boson, the 
genuine two-loop calculation is particularly simple in that gauge. It is, however, instructive 
to inspect the situation in general gauge for the Z boson as well, a task that is pursued 
in the Appendix. 

The reducible two- loop graphs involve the product of the one-loop induced gg*H and gg*Z 
vertices and can be calculated with conventional one- loop techniques. These graphs by 
themselves form a gauge-invariant, UV- finite and IR-finite subset of diagrams. 

In practice, we have performed two completely independent calculations of all loop contri- 
butions. Version 1 follows basically the strategy described in Ref. [16] for the related pro- 
cess of scalar-pseudoscalar Higgs-boson pair production. Here the Feynman diagrams are 
generated with the program FeynArts [33] , and the large-mass expansion of the diagrams 
involving top-quark loops is performed by inhouse Mathematica routines. The calcu- 
lation is carried out in Landau gauge, so that no two-loop diagrams with massless-quark 
loops contribute. In the second calculation, the diagrams are generated by QGRAF |34| 



and expanded using EXP / Q2E 35 , 36 . The massive two- loop vacuum integrals resulting 
from top-quark loops are calculated using MATAD [37j. This calculation is carried out 
both in Landau and in unitarity gauge. In the latter, the Goldstone bosons are absent, 
but massless-quark loops contribute to the ggZ* vertex; these graphs are calculated with 
the program MINT |38|. While the results of the two different calculations in Landau 
gauge are in mutual agreement term by term, for the calculation in unitarity gauge only 
the full result agrees (after the renomalization procedure). 

Including all required counterterms, the virtual contribution is given b}Q 

Vvirt = J dPS 2 [|A^ (e)| 2 (l + 25Z + 5 CS ) + 2Re \ji^e)MoW}] , (5) 

where A4o(e) is the LO amplitude in D = 4 — 2e dimensions, Ai\(e) the amplitude of 
the virtual corrections at NLO, and (IPS2 denotes the 2-particle phase-space element. We 
renormalize the strong coupling in the MS scheme assuming nj = 6 flavours, and the 
top-quark mass in the on-shell scheme. Note that the renormalization factor is gauge 



1 Note that we consistently suppress terms involving ln47r and je, which accompany poles in e, because 
they cancel in the UV-finite result as usual. 
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dependent; in Landau gauge, it is 



8Z = ZlZzZJ^Zr. 



25Z g + 5Z 3 



2e5Z m + SZF 



5Z„ 



6Z* 



5Zr, 



6Z? 



TT 



5 ! 

1 

4? 



a Sm 1/1 , 
-_iT R - -+ln 
tt 3 \ e 



m: 



(6) 



TT 



c F - + o(r 



-2— C F , 

TT 



where Cf 



4 t 



i, and Ca = 3 are the QCD colour factors. In unitarity gauge, 



the mass counterterm 5Z m is absent, but 5Z§ = — ^Cf |32| is needed for the ggZ vertex 
instead of 6Z$. 

The term <5cs, 



<5cs 



as //r 

TT V S 



C A 



+ 



67 
18 



TT 

~6 



+ 



C A 



+ 1 



11 



10 



Trui 



-C a - -Trui 



TT 

( 

12 



(7) 



in Eq. §5§ needs to be added according to the dipole subtraction method [39] , where ni = 5 
denotes the number of light flavours. 

In order to be consistent with the currently available PDF sets, we express our final result in 

(5) 

terms of a s , the strong coupling with five active flavours using the matching relation 



40 



a. 



(5) 



(/*) 



tt. 



(6) 



(aO 



i + 



tt : 



(6) 



(m) 



7T 



1 



In 



+ 0(a s 2 



(8) 



so that the logarithmic dependence on m t vanishes. This procedure is equivalent to de- 
coupling the top quark in the a s renormalization by subtraction of the top-quark loop in 

the gluon self energy at zero-momentum transfer, instead of using the MS prescription. In 

(5) 

the remainder of this paper, we set a s = a s . 

Inserting the QCD colour factors, the final result for the virtual contribution can be written 

as 



°Virt 



1 + 



Wf 164 23 
V 9 6 



TT 



Clo + <?"(virt,red) j 



where 



<tlo = J dPS 2 \M \ 2 , 



(9) 



(10) 
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(cf.Eq.(gJ() and 

°~ (virt,rcd) denotes the contribution from the type of reducible diagrams 
shown in Fig.[2](k). It is not proportional to cr LO ; as a function of the partonic Mandelstam 
variables s = (pi + P2) 2 , t = (pi — pz) 2 , and u = [pi — £>h) 2 > it reads 



a s \ 3 a 2 7r 2 I 



^(virt.red) " / ^2 ^ - J 76gs 4 £ 4 M | g _ M | 



x < (sM|m£ + sM| - s 2 M£ - 2s 2 Mf + s 3 ) 



/-u\ M% , f-i\ M% \ (11) 
2 + ln[ — A- + In f — ^ 1 1 



VMj/u-Mf \Mpi-Ml 
+ ( - sM|m£ + sM| + s 2 M^ - 2s 2 M| + s 3 ) 

^ H ^ + In — ^ 7 o\~<t + In I — - l —. ^ 



i-Ml u-Ml \Mp {u-Mlf \Ml)(i-MlY 
(in) NLO real corrections 

The real corrections are induced by the partonic channels gg — > HZg, gq — > HZg, gq — > 
HZg, and qq — > HZg, where in the channels involving external quarks only the squares 
of the diagrams with closed quark loops are taken into account. At first sight, the most 
complicated one-loop diagrams are pentagon graphs with a top-quark loop. However, 
in the large-top- mass limit, these graphs vanish. The algebraically most complicated 
diagrams are the box graphs with external Z*ggg fields; they are the only ones that receive 
contributions from the vector-coupling of the Z boson, while all other diagrams (summed 
in pairs of opposite charge flow in the loop) depend only on the Z-boson axial-vector 
coupling. 

The actual calculation of the diagrams can be performed using standard one-loop calcu- 
lational techniques and has been carried out in three completely independent ways. The 
first approach builds on graphs from FeynArts 1.0 1 33 1 and reduces or expands the full 
amplitudes with inhouse Mathematica routines, which produce output in the form of 
Fortran code. The occurring one- loop tensor and scalar integrals are numerically evalu- 
ated with the (not yet public) library Collier that is based on the results of Refs. [41 -46 



The second calculation is based on the program packages FeynArts 3.2 47 and Form- 
Calc/LoopTools |48,49], as far as the calculation of one-loop graphs is concerned that 
do not involve top-quark loops. The large-mass expansion of the top-quark loops here 
again is carried out using inhouse routines (independent from the ones of version 1). 

The third approach again uses the QGRAF/EXP/Q2E/MATAD-setup for the generation and 
expansion of the diagrams and the evaluation of the massive vacuum diagrams. The cal- 
culation of the massless triangle and box diagrams is performed by an extended version of 
the FORM 1 50 1 routine previously used in Ref . [ll] , which implements algebraic Passarino- 
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Veltman reduction [42] and analytic results of the scalar integrals given in Ref. |51| . Here 
we explicitly verified the gauge invariance of the result with respect to the Z propagator 
as well as to the external gluons. For the latter, we assumed a general axial gauge, where 
the polarization sum reads 

Effiforfffo) = *"» [-sr + pl X + f n " ) (i = 1,2,3) (12) 

with an arbitrary light-like vector n which drops out in the squared amplitude. 

All real-emission channels contain IR singularities in their integration over phase space. 
More precisely, gg —> HZg becomes IR singular if the emitted gluon becomes soft or 
collinear to one of the incoming gluons. The other channels involve only collinear singu- 
larities. The separation of the IR singularities in the phase-space integration is achieved 
using the standard dipole subtraction method [39 1, where an auxiliary cross section is 
subtracted from the full real-emission part and added back after an analytical integration 
(in D dimensions) over the one-particle emission phase that contains the IR singularity 
(cf.Eq.0). 



4 Numerical results 

4.1 Input values 

We use the following input parameters: 

M z = 91.1876 GeV, M w = 80.399 GeV, = 1.16637 • 1(T 5 GeV~ 2 , 
m b = 4.75 GeV, m t = 172GeV. 

In the effective-theory approximation, we set m t — > oo and rrtb = 0, as discussed above. 
For the electromagnetic coupling constant a we employ the scheme, where the coupling 
constant is defined as 

«- ^ , S w -L C w -1 (14j 

As the default PDF sets, we use MSTW2008(N)LO [52] when evaluating a (N)LO quantity. 
The corresponding input values for the strong coupling are given by a\°{M'i) = 0.13939 
(a^ LO (Mz) = 0.12018). The running of a s is performed to the order under consideration: 

^ = -i>s +2 A, as) 

with n = (n = 1) at LO (NLO). Both the PDFs and the a s evolution are implemented 
with the help of the LHAPDF library [53]. 
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Figure 3: Comparison of the LO partonic cross section in the effective (labelled "LME") 
and the full theory for various, partly hypothetical values of the top-quark mass. The 
curve for mt = 3.44 TeV cannot be distinguished from the LME result. 



Our default choice for the renormalization and the factorization scales /ur and hf is the 
invariant mass of the HZ system: 

IM) = Vipu+Pl) 2 ■ (16) 



4.2 Leading-order considerations 

In order to get an idea about the quality of the effective theory, we show some studies 
at LO before presenting our NLO results. Figure[3] shows the partonic cross section both 
for the exact top-mass dependence and in the effective theory. The exact result exhibits 
a kink at the top-quark pair threshold = 2mt = 344 GeV which clearly cannot be 
reproduced by the effective-theory approach. For larger values of we do not expect an 
expansion in 1/mt to converge. In fact, higher-order terms in this expansion would most 
likely worsen the prediction in the region of larger yS. 
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120 140 160 

M H [GeV] 
(a) Inclusive cross section 



M H [GeV] 
(b) p T ,H > 200 GeV 



Figure 4: Comparison of the LO hadronic cross section in the effective and the full theory 
for ^/s = 8TeV (dashed) and 14TeV (solid). 



Taking into account the kinematical constraint Vs > Mh + M%, the region where the 
effective theory is nominally applicable shrinks to zero for Mh > 2m t — M% ~ 253 GeV. 
Figure[4](a) compares the total inclusive LO hadronic cross section at 8TeV and 14TeV 
when the full top- and bottom-mass dependence is taken into account to the effective- 
theory result. The behaviour is expected from the considerations above: The effective 
theory works better for smaller Higgs masses, agreeing to the full results within 2% (25%) 
for 8TeV (14TeV) at M H = 125 GeV. Note that the PDFs suppress the contribution from 
larger s, thus emphasising the region where the 1/mt expansion converges. For larger 
values of Mh, the effective-theory approximation deteriorates; at Mh = 200 GeV, the 
deviation to the full result is 74% (143%) for 8TeV (14TeV). 

The situation becomes more problematic in the boosted regime which we study by im- 
posing a lower cut on the Higgs' transverse momentum, requiring pr,H > 200 GeV, see 
Fig.[|(b). In this case, the minimal value for Vs is already above the top-quark thresh- 
old when Mh = 100 GeV. Consequently, the direct application of the effective-theory 
approximation is off by almost a factor of five to ten, which is clearly unacceptable. 

A direct evaluation of the NLO contribution in the effective theory is therefore not possible. 
However, in Refs. 21-27 it was shown for the process gg — > H at NLO and NNLO that 
the perturbative correction factor, defined at NLO in Eq. Q, depends only very weakly on 
the top-quark mass. To some degree, this holds even far outside the convergence region of 
the heavy-top expansion, as long as only the leading term in 1/mt is taken into account. 
Motivated by this observation, we move on to NLO and present our results in the next 
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Figure 5: NLO hadronic cross section as obtained by using Eq. (17) (upper), and NLO 
K-factor (lower) for yfi = 8TeV (dashed) and 14TeV (solid). 



section. 

4.3 Next-to-leading order results 
4.3.1 Correction factor 



As outlined above, we evaluate the NLO hadronic cross section by rescaling the full LO 
result by the perturbative JC-factor calculated in the effective theory: 

CT ap L p ?ox( m t,m b ) = cr LO (m t ,m h )K(m t -> oo,m b = 0) 

a LO (m t ,m h ) NLO (17) 

a (m t — > oo, mb = 0) . 



cr LO (m t — > oo,mb = 0) 
Since we are aiming at a NLO quantity, it actually might be more appropriate to evaluate 



the formally LO cross sections in Eq. (17) with NLO PDFs. We checked that the effect 
of this is much smaller than the uncertainty due to variations of the renormalization and 
factorization scale, which is why we stick to LO PDFs in a h °. 
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Figure 6: Individual contributions to the NLO hadronic K-factor at \/s = 14TeV as 
discribed in the main text. 



Figure [5] shows the gluon-induced cross section obtained in this way for sfs = 8 TeV and 
y/s = 14 TeV hadronic centre-of-mass energy, together with the corresponding perturba- 
tive correction factor K. Part (a) of Fig. [5] shows the total inclusive cross section, while 
in part (b) the boosted scenario with pt,h > 200 GeV is shown. In both cases, we ob- 
serve a K-factor of the order of two, almost independent of Me, with a slight increase 
towards lower centre-of-mass energies. This behaviour is very similar to the one observed 



for gluon-induced single-Higgs 21 - 27 and Higgs pair production |16|. The correction even 
slightly exceeds the well-known correction factor for gg — > H. 

A breakdown into individual contributions to K is shown in Fig. [6] for the total inclusive 
cross section at \fs = 14 TeV: 

• i^LO — change of PDF sets from LO to NLO 

• Aifggvirt — virtual corrections including integrated dipole terms according to Ref. [39] 

• A.fr ggrea i — correspondingly regularized real corrections 
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0.1 



(a) Inclusive cross section 



M/MO 
(b) p T ,H > 200 GeV 



Figure 7: Scale dependence of the hadronic LO and NLO cross section for y/s = 8 TeV 
(dashed) and 14 TeV (solid). The renormalization and factorization scales are varied 
simutaneously around the central scale fj,Q = \J (pn + Pz) 2 - The Higgs mass is set to 
M H = 125 GeV. 



• Aif q g, AK m — contributions from qg and qq initial states 



The sum of all these terms results in Ktot, the total K- factor. 



4.3.2 Residual scale uncertainty 

As described in the introduction, the LO scale dependence for this purely gluon-induced 
process is quite large. NLO corrections typically decrease this uncertainty. Let us recall 
the situation in the gluon-fusion process gg — > H, however: For the LO result, the usually 
adopted scale variation by a factor of two around the central scale leads to a gross un- 
derestimation of the size of the higher-order effects. At NLO, the scale uncertainty is not 
significantly smaller than at LO, but it does provide a good estimate of the NNLO effects. 
Consistently, inclusion of the NNLO corrections leads to a significant reduction of the scale 
uncertainty. 

Expecting a similar behaviour for the gg —> HZ process, it is not surprising to see the result 
shown in Fig.[7j Both for the inclusive and the boosted scenario the scale dependence 
decreases from more than 100% at LO to 60% at NLO when the renormalization and 
factorization scales are varied simultaneously by a factor of six around their central value 



Ho, see Eq. (16). As for the process gg — > H, the behaviour in h/hq is strictly monotonous, 
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and the LO and NLO curves do not intersect. Therefore, a preferred value for and ^r 
cannot be deduced from these plots. The radiative corrections increase with ///MO) so there 
is a slight tendency towards choosing smaller values of \i. Nevertheless, in our numerical 
analysis we stick to the "natural" value \x§ as the central choice. 

Note that, also similar to what is observed in gg — >■ H, variation by a factor of two would 
not lead to any overlap between the LO and the NLO predictions. 

The similarity between the processes gg — > H and gg — > HZ suggests that the NLO error 
estimate due to scale variation is quite reliable for the process gg —> HZ. In order to take 
into account the fact that the effective theory is expected to work not quite as well in 
gg — > HZ as in gg — > H, we determine this uncertainty by varying /i within a factor of 
three rather than two around the central value ^o- The numerical results are listed in 
TableE 



4.3.3 Total inclusive cross section 

In this section we provide the most up-to-date numbers for the total inclusive cross section 
for the Higgs-strahlung process at the LHC with 8 and 14TeV, including 



• NNLO Drell-Yan terms a HV ' DY of order 5 4 a™ (n = 0, 1, 2) [?}]9]; 

• electroweak corrections which are applied as an overall factor to the Drell-Yan terms 



12,14 



• top-loop-induced corrections of order 0(Xtg 3 ag) 

• gluon-induced terms of order A^<? 2 a" (n = 2,3); n = 3 corresponds to the newly 
calculated terms of this paper. 

For the non-gluon-fusion part of the cross section, the scale variation is obtained by using 
the MSTW2008NNLO PDF set and varying /ip and ^r independently within the interval 



,(-) 



(+) 



(a*r,/xf)/M h € [1/3,3] x [1/3,3], which results in a cross section interval [<7 no _gg, cr no . 
The central value of the total cross section is then obtained as 



^central 



a (+) +a (+) +a (-) +a (~) 
^ no-gg ' gg 1 no-gg ' gg 



(18) 



where <7gg ^ are the boundaries of the scale uncertainty interval of the gluon-induced compo- 
nent which can be obtained from Table[TJ Accordingly, the scale uncertainty is calculated 



as 



A 



scale 



a (+) +ff (+)_ (7 (-) 

>-* j^O-gg ' gg 



no-gg "gg 



/(2a, 



central I 



(19) 
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fs [TeV] M H [GeV] 



a, 



NLO 

If 



[fb] 



no Pt,h cut 



8 


115 


8 


120 


8 


125 


8 


130 


14 


115 


14 


120 


14 


125 


14 


130 



p T ,n > 200 GeV 



19.8 
18.7 
17.7 
16.7 



+61% 
34% 

+61% 
34% 

+61% 
-34% 

+61% 
34% 



70 1+51% 
' a -31% 

+51% 
-31% 

+51% 
31% 

+51% 
31% 



75.1 
71.1 
67.2 



39.3 
37.2 
35.1 
33.1 
152 



144 
136 
129 



+32% 
24% 

+32% 
24% 

+32% 
-24% 

+32% 
-24% 

+27% 
-21% 
+27% 
-21% 

+27% 
-21% 

+27% 
-21% 



8 


115 


1 41+65% 
1.41_ 36% 


2 94+34% 
Z - y ^-25% 


8 


120 


i oq+65% 
1.33_ 36% 


9 70+33% 
Z -' y -26% 


8 


125 


i 9R +65% 
±.zu_ 36% 


9 RQ+34% 


8 


130 


1 iq+65% 
i.iy_ 36% 


2.48+33% 


14 


115 




i q o+29% 


14 


120 


a r:q+55% 


i o i +28% 

io.i_ 22% 


14 


125 




19 c+29% 


14 


130 


5.87+^ 


11 S+29% 

ii.o_ 22% 



Table 1: Cross sections of HZ production via gluon fusion for LHC energies in the range of 
phenomenologically preferred Mh values. The scale uncertainty is given in percent. The 
latter results from a rescaling of /xr = [ip by factors of 3 and 1/3 relative to /j,q. 
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y/s [TeV] 


M H [GeV] 


<r[pb] 


A SC alc[%] 


ApDF+ct s [%] 


A to tal[%] 


pp -> HW 


8 


115 


0.926 


±0.6 


±2.3 


±2.9 


8 


120 


0.805 


±0.6 


±2.5 


±3.1 


8 


125 


0.705 


±0.6 


±2.3 


±3.0 


8 


130 


0.617 


±0.7 


±2.4 


±3.1 


14 


115 


1.97 


±0.6 


±2.0 


±2.6 


14 


120 


1.73 


±0.7 


±1.8 


±2.5 


14 


125 


1.52 


±0.7 


±2.2 


±2.9 


14 


130 


1.34 


±0.6 


±2.0 


±2.6 


pp -»■ HZ 


8 


115 


0.540 


±2.6 


±2.4 


±5.0 


8 


120 


0.475 


±2.8 


±2.4 


±5.1 


8 


125 


0.419 


±2.9 


±2.4 


±5.3 


8 


130 


0.371 


±3.1 


±2.3 


±5.4 


14 


115 


1.24 


±3.5 


±1.8 


±5.3 


14 


120 


1.10 


±3.7 


±1.6 


±5.3 


14 


125 


0.983 


±3.9 


±1.6 


±5.4 


14 


130 


0.880 


±4.0 


±1.9 


±5.9 



Table 2: Total inclusive cross section for the processes pp — > HW and pp —> HZ. The 
latter includes the newly calculated NLO gluon-induced terms. The evaluation of the scale 
and PDF±a s uncertainties is described in the main text. 



The influence of the newly evaluated NLO gluon-induced terms on the overall PDF±a s 
uncertainty is rather small, since this contribution comprises only about 5% of the total 
cross section, and will be neglected. Therefore, we base the estimate of the PDF±a s 
uncertainty solely on what is currently contained in vhOnnlo (i.e., LO gluon-induced terms 
are taken into account). Following the PDF4LHC recommendations |54| by using the NNLO 
PDF sets from MSTW2008 [52], CT10 [55], and NNPDF23 [56], we obtain a cross section 
interval [c PDF+as , o"pdf+oJ an< ^ calculate the resulting uncertainty as 

A _ °PDF+a s °PDF+a s / on \ 

A PDF+as - — ^ • (20) 

°PDF+a s ' °PDF+a s 

Our results are shown in Table[2j We find that the NLO gluon-induced terms calculated 
in this paper increase the central values of the HZ cross section by about 4% (7%) at 
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8TeV (14TeV). Since the K-factor for these terms is of the order of two, and their scale 
uncertainty decreases by roughly the same factor when going from LO to NLO (see TableHJ, 
the overall scale uncertainty on the total inclusive cross section remains almost unaffected 
by the inclusion of the new terms. For completeness, we also include updated numbers 
for HW production in Table[2j even though they are not affected by the NLO gluon-fusion 
terms calculated in this paper. 

As a side remark, we note that the PDF+a s uncertainties of Table[2] are significantly 
smaller than in Ref. [11]. This is due to the use of only NNLO PDF sets in this newer 
version, while the previous numbers were based on a rescaling of the NNLO MSTW2008 
uncertainty by the NLO PDF error. For HW production, also the scale uncertainty is 
slightly smaller in Table[2]than in Ref. |11| . This is because these previous numbers were 
obtained by linearly adding uncertainties of the "top-induced" terms to the rest, while 
here we vary the scale in both contributions simultaneously. For the HZ cross section, this 
is overcompensated by the uncertainty of the NLO gluon-fusion component, see above. 

5 Conclusions 

The gluon-induced corrections to the Higgs-strahlung process have been calculated through 
NLO, i.e. O(ag). The perturbative correction factor is obtained in the limit nit — > oo and 
mb = 0, and then used to rescale the full LO cross section in order to obtain a prediction 
for the NLO result. The behaviour of the perturbative series and the residual scale vari- 
ation are found to be comparable to the gluon-fusion processes of single-Higgs and Higgs 
pair production. The success of the effective-field-theory approach in describing correc- 
tions to single-Higgs productions, where exact calculations for higher-order corrections 
are available, gives confidence in a reliable prediction of the theoretical uncertainty due to 
higher-order effects for gluon-induced HZ production. 

Numerical results were provided for current and future LHC energies, both for the fully 
inclusive cross section and for boosted Higgs kinematics. We use these results in order 
to provide the most up-to-date predictions for the hadronic Higgs-strahlung process at 
relevant collider energies. The large corrections on the gluon-induced terms, combined 
with their large scale uncertainty, increases the overall uncertainty on the total HZ cross 
section by about 1%. 

In the near future, the results will be included in the publically available numerical program 
vhOnnlo. 
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Appendix 



Treating gg — > HZ via the anomaly relation 



In this appendix we sketch yet another way to calculate the genuine two-loop part upon 
employing various tricks, such as the Landau- Yang theorem 17 [18 , the Adler-Bell-Jackiw 



(ABJ) anomaly relation 57 , and the Adler-Bardeen theorem [58 , in order to reduce the 
calculation to a much simpler massless one-loop calculation. For definiteness we employ 
the conventions of Ref. [59] for the SM parameters and Feynman rules. 



In Section[3]we explained that the LO and NLO virtual amplitudes for the process gg 
have the following properties: 



HZ 



• Non- vanishing contributions originate only from loops of vertex types ggZ, ggG , 
and ggH, as illustrated in Fig.[2|L,c,k,l,m), while all graphs with more than three 
external legs attached to the loop vanish or compensate each other for wit —> oo and 
m qj L t = 0. 

• At NLO the loop-induced ggH vertex is only relevant at the one-loop level within 
reducible diagrams like the one shown in Fig.[2^k). Owing to their simplicity we do 
not consider those reducible diagrams in the following. 

• When the ggZ vertex is attached to an intermediate Z propagator, and only in this 
case two- loop diagrams of this vertex become relevant, the Landau- Yang theorem 
implies that only the longitudinal part of the Z propagator contributes. 

Thus, the contribution of one-particle-irreducible (1PI) diagrams, i.e. the graphs in Fig. [2] 
(b,c,l,m), to the matrix elements Ai n (n = 0, 1) is given by 

M]f\e) = Gfc&fan, -k) Df(k) ^ e\ a 

s w c w (21) 

- ie?4 Gf/ G °( Pl ,p 2 , -k) D^k) (pz + 2 m )° ^ 

where we use an obvious notation for the external polarization vectors £i,£2, £g for the two 
incoming gluons and the outgoing Z boson, a, b are the gluonic colour indices, k = P1+P2 is 
the momentum of the intermediate Z or G lines, G^Jp [pi,P2, —k) and Gfn, [Pi,P2, —k) 
are the amputated Green functions for the ggZ/G° vertices, and D^ a (k) and D^(k) are 
the Z/G° propagators in the general gauge, 

D f< fc >= zi fef -'^rar D ^ ) = rar (22) 



22 



Inserting the latter into Eq. (21 ) and exploiting the fact that the transversal part of D ( T{k) 



(first term in Eq. (22)) does not contribute, we obtain 

e Oh • 4) 



l c 2 



i^f ( M , -k) ^ - Gf/ G °( PuP2 , -k) 



(23) 



where we have used k = pz + pn, Pz ■ £z = 0, and /c 2 = s. The terms in square brackets 
can be simplified using the well-known form of the ABJ anomaly. 

To this end, we recall that the ABJ anomaly relation for the axial current = ip fj^jst/jf 
and the pseudo-scalar operator p/5 = ipf^ipf for a quark q reads 



d »3q 5 = 2[m qPq,5 



n. 



4tt ^ ' 



(24) 



where F a '^ = \^ vp(J F^ a is the dual of the gluonic field-strength tensor F a '^. Eq. pi} 



is an operator relation, valid for bare quantities, expressing chiral symmetry. The Adler- 



Bardeen theorem 58 states that it is correct to all orders in regularization schemes that 



respect chiral symmetry. In regularizations that are not chirally symmetric, Eq. (24) has to 



be restored by extra counterterms from evanescent operators. For the 't Hooft-Veltman 75 
scheme 30i 31 , which is employed in our calculation, these are the counterterms 5Z§ and 



5Z^ calculated in Ref. [32] and already used in Section[3l The relevant momentum-space 
Feynman rules for the composite operators j^ 5 , p q ^, and FF are given by 




The dotted line in the Feynman rules indicate that the momentum p flows into the vertex. 
The operator FF induce Feynman rules with three and four gluon lines as well, but those 
will not contribute in the following. The fermionic Feynman rules are related to the 
couplings of Z/G° to the quark: 

gj3 

Z»qq ■ -^^7^75 + vector part, G ^ : ^75, (26) 

zs w c w s w 

where is the third component of the weak isospin of q. 

Now we can make contact with the Green functions G g ° gi>z and G gag6G ° introduced above. 
Since we consider only QCD corrections, in the relevant graphs the couplings of the external 
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Z/G° lines to an internal quark line represent the only electroweak coupling in the ggZ/G° 
vertex functions. These electroweak couplings can be interpreted as the insertions of the 
operators jf* 5 and p q ,5, because the vector part of the Z coupling does not contribute, as 



explained in Section[2j Thus, we obtain 

i 

r g a g b c° sr^ r g a g b 



b G c 



ng a g b 1 



rs a g b G° 



le/ w,g n g a g b j q , 5 



<-J;j 



(27) 



s w c w M z 



which is valid up to NLO QCD, where G gagb7j ^ G ° \ q denotes the contributions induced by 
closed loops with quark q. Using this relation, Eq. (24) implies: 



\p H 



-M z Gl a / G °( Pl ,p 2 ,p) 



+ 



irT "- <¥ ' a G&WfaM), (28) 



2s w c w 4tt ^ 



where we have restored the arguments of the incoming momenta. Keeping in mind that 
p = —k, we thus can write the contribution of closed q-loops to M.^ l {e) in Eq. (23) in two 
different ways: 



M 



XPli 



e (PH • 4) Ji 



e7e 



1°2 



e Oh • 4) 



ie/. 



+ 



<z 2s w c w 4vr s-£M| 



■ e7e 



1°2 



(29) 



q + 2s w c w 4m s — £,M£ 



Now we exploit our approximation of mt — > oo and m^ t = 0, which implies k p G gag z |t = 
and G gagbG °\ qj tt = 0. The former result is taken from our diagrammatical large-mass 
expansion, the latter is a trivial consequence of the vanishing Yukawa couplings of the 
massless quarks. For the top quark we take the second form of the last equation of 
Eq. (29) and for the other quarks the first, and obtain 

a 



g (PH • g|) u c v 

s w c w s 

ie -^w,t a s s 



2s w c w 4vr s - £M| iPl,P2 ' } 



+ 



2 GlT FF \ Pl ,P2,-k) 



2s w c w An s — £M Z 

e^ 2 Gf/( FF \ Pl ,p 2 ,-k) 



(30) 
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Figure 8: NLO QCD diagrams for the Green function G gagb ( FF \ 



where the dependence on the gauge parameter £ cancels, as it should be. For the LO 
matrix element, the Green function G g<lg ( FF ^ just has to be replaced by the Feynman rule 
for the FF operator with two external gluon legs, yielding 



Mr 



-5 



(PH • 4) 



S W C W 



M z 



(31) 



in agreement with Eq. ([3]). The NLO QCD corrections to G gag ( FF ") are induced by the 
diagrams shown in Fig.[8j The actual calculation of these one-loop diagrams is very simple. 
For on-shell gluons only the first two diagrams contribute and yield 

' ,0 " ~ - sC (s, 0,0, 0,0, 0,0) 



lPI.virt 



(0 



2tt 



(32) 



where Co is the one- loop scalar 3-point integral in the convention of Refs. 44 -46] (see e.g. 
Ref. 1 60 1 for the explicit result). This is the unrenormalized result for the virtual correction 
which still has to be renormalized. As stated above, the anomaly equation (24) is valid 
for bare quantities, so that the FF term receives renormalization contributions from a s 
and the gluon field. Since in our case only the two-gluon contribution of FF is relevant, 
we get the simple factorizing contribution to the NLO amplitude, 

M\ Pl ' ct (e) = (25Z 9 + 5Z 3 )M . 



(33) 

Combining the renormalized virtual amplitude _A/fJ PI,virt +7W} PI,ct with the contribution 
5qs j see Eq. Q , of the Catani-Seymour /-operator of the subtraction function yields the 
1PI part of the correction to the cross section, as given in Eq. ([9| by first term on the r.h.s. 
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